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Abstract 

As a metric for amplitude fluctuation of orthogonal frequency division multiplexing (OFDM) signal, 
cubic metric (CM) has received an increasing attention because it is more closely related to the distortion 
induced by nonlinear devices than the well-known peak-to-average power ratio (PAPR). In this paper, 
the properties of CM of OFDM signal is investigated. First, asymptotic distribution of CM is derived. 
Second, it is verified that 1.7 times oversampling rate is good enough to capture the CM of continuous 
OFDM signals in terms of mean square error, which is also practically meaningful because the fast 
Fourier transform size is typically 1.7 times larger than the nominal bandwidth in the long-term evolution 
(LTE) of cellular communication systems. 


Index Terms 


Cubic metric (CM), fast Fourier transform (FFT), orthogonal frequency division multiplexing (OFDM), 
oversampling, peak-to-average power ratio (PAPR). 


I. Introduction 

Orthogonal frequency division multiplexing (OFDM) is an attractive multicarrier modulation 
technique for broadband wireless access systems due to its strong immunity to multipath fading 
and high spectral efficiency. However, OFDM signals suffer from high amplitude fluctuation 
which causes performance degradation due to nonlinear devices. A well-known metric for 
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amplitude fluctuation of OFDM signal is peak-to-average power ratio (PAPR). Many research 
efforts have been carried out to find efficient PAPR reduction techniques [[Il-[[6l. Also, the 
distribution of the PAPR of continuous OFDM signal was derived [|7]| and it is widely accepted 
that four times oversampling is enough to capture the PAPR of continuous OFDM signals [[H. 

Another metric for amplitude fluctuation of OFDM signals has been considered ll9lffT0l . which 
is known as cubic metric (CM) ffTTl . Studies on PAPR and CM suggest that, except for large 
power backoff, CM is more closely related to the amount of distortion induced by a nonlinear 
power amplifier than PAPR [fTOl . Moreover, after analyzing certain OFDM-type signals that are 
considered to meet the goal of the long-term evolution (LTE), it was shown in the 3GPP that 
CM predicts amplifier power de-rating more accurately than PAPR [fT^ . 

Thus, recent research to reduce the CM for the LTE systems has been carried out [|T3ll - [|T^ . 
For example, in ffT3ll . the clipping and filtering method to reduce the CM instead of PAPR is 
proposed, where the descent clipper different to the conventional clipper is designed. In [[T4l . 
the tone reservation method is proposed, where minimizing the CM value is set to the objective 
function of convex optimization. 

A great deal of the literature has been devoted to research on the properties of the PAPR 
metric as in [|7l|[l8l|. On the contrary, such analysis on the CM has not been done before. Thus, it 
is worth revealing more about the behavior of this metric. In this paper, an asymptotic probability 
distribution of CM for continuous OFDM signals is derived. Also, sufficient oversampling rate 
for capturing the CM of continuous OFDM signals is obtained. 


II. Cubic Metric 

The CM of OFDM signals is defined as IfT^ 

RCMl 


CMLe = 


IdB 


RCM 


ref 


dB 


K 


where RCM is the raw CM defined by 


RCM[s(f)]|^g^201og 


rms 


Ht)\ 


( 1 ) 


for a continuous OFDM signal s{t) and both RCMref 
considered OFDM systems flT^ . As an example, in the downlink of LTE, RCMref 


\rms[s(f)] 

and K are determined according to the 

= 1.52 dB 

dB 
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and K 


1.56 are used. Thus, from ([T]), we are only interested in 


RCM[s(t)] 



where Pav is the average power of the eontinuous OFDM signal s{t). To simplify analysis, we 
will use the square of it as 


^ = (RCM[s(f)])2 = E 


m\ 

V Pav 




In praetiee, instead of ealeulating the CM of eontinuous OFDM signals, we ealeulate the CM 
of diserete OFDM signals. Let s{t) be a eontinuous OFDM signal and its L times oversampled 
OFDM signal sequenee be represented as 


Sn,L — s{nTs/LN), 0 < n < LN — 1 


where is the OFDM signal period, N is the number of subearriers, and L is a real number 
larger than or equal to one. 

Without loss of generality, the input symbols in frequeney domain are assumed to be statis- 
tieally independent, identieally distributed (i.i.d.) random variables with zero mean, where the 
input symbol is the eomplex data of eaeh subearrier. Then the OFDM signal eomponents in time 
domain are given by the sum of i.i.d. random variables. Thus, from the eentral limit theorem 
(CLT), the magnitude of Sn,L is Rayleigh distributed [|3- Therefore, if it is normalized as 

r{t) = \s{t)\/^/P^ 

'^n,L — I'StIjLI/ \/Pav) 

the probability distribution funetions (PDFs) of r{t) and rn,L are given as 

fr(t){r) = = 2re-^\ 

Finally, for the diserete OFDM signal sequenee obtained by L times oversampling, ^ and RCM 
are expressed as 

^ LN-l 

^ Tn ^ 

n=0 

and 

RCMl = y/Ti- 
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III. Properties oe r® ^ 

In (|^, is sample mean of r® ^’s and thus in this seetion we investigate the properties of 
the random variable By the definition of Weibull distribution, the power transformation 
Wn,L — ^nL Rayleigh distributed random variable rn,L is known as Weibull distribution 

uni. The PDF and eumulative distribution funetion (CDF) of Wn,L are given as 

= ^w“^exp(-w;^) 


and 

= 1 -exp(-'«;5), 


respeetively. The /cth-order moment of Wn,L is known as = r(l + 3k), where r(a) = 

(a — 1)! is the Gamma funetion for an integer x and E[-] denotes expeetation value. Then we 
have 

EK,l] = 3! = 6 (3) 

EkiJ = 6! = 720 

VaT{wn,L) = 684 (4) 


where Var(-) denotes varianee. 

Now, we ealeulate the eovarianee of two random variables Wn,L and Wn',L', whieh will 
be denoted as CoY{wn,L-,Wn',L')- For this, we first obtain it with the eontinuous time lag r, 
Cov{w{t),w{t + r)), where w{t) = r^{t). The joint moment of w{t) and w{t + r) is expressed 
as uni 


+ r)] = (1 - p2)i+3p+3<?r(i + 3 p) + 3 ^) + 3p,l + 3q; 1; pi) (5) 


where 2 -Fi(-,-;S') is the Gauss hypergeometrie funetion and Pr = Ps{t),s(t+T) is the Pearson’s 
eorrelation eoeffieient between s{t) and s{t + r) ffTTl . From ([^ and (|^, we have 

Cov{w{t),w{t + t)) = Elw(t)w(t + r)] — E[w(f)]E[M;(f + r)] = 36(9p^ + 9p^ + pi). 

It ean be assumed as in 0 that the power speetrum of baseband OFDM signal has eonjugate 
symmetry, whieh is valid beeause the power speetrum of baseband OFDM signal ean be designed 
to have symmetry at the eenter of the bandwidth by giving proper frequeney offset. This 
assumption guarantees the autoeorrelation funetion of s(t), E[s(f)s(f + r)*], to be a real funetion. 
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In any case, the frequency offset is immaterial to our analysis, beeause it does not ehange the 
magnitude of the envelope of OFDM signal. In this ease, the normalized autoeorrelation funetion 
Pr is known as m 

sm{7rNT/Ts) 

N sin(7rr/Ts) 

for r 7 ^ 0 and elearly p,- = 1 for r = 0. Note that Pr = 0 when 

T, 2T, 3T, 

r = ±— ±—^ ±—^ ■ ■ ■ 

AT’ iV ’ iV ’ 

whieh implies the well known faet that the elements of the Nyquist sampled OFDM signal 
sequence are mutually independent. 

The time lag between two diserete samples Wn,L and Wn',L' is 

riTg n'Tg 

LN~ Un' 


Finally, the eovarianee of Wn^L and tCn'.L' is given as 


Go\{Wn,L,Wn',L') 


36(9p2+9p4+^6) 


n'Ts 

L'N 


( 6 ) 


IV. Distribution of RCM 

In this seetion, we obtain the asymptotie distribution of RCM^, by investigating the distribution 
of first. 


A. Mean and Variance of 

The mean of is elearly 6 from Q beeause is sample mean of Wn,L- To find the varianee 
of suppose that LN is an odd integer. Even though LN can be any real number larger than 
or equal to N, it is not diffleult to show that diserepaney by the assumption is negligible. Sinee 
Sn,L is a eomplex stationary Gaussian proeess, both and Wn,L are also stationary random 
process. Therefore, the varianee of beeomes 

Var(^i) =al = ^ ^ Cov(w;o,l, wpi). (7) 

We can separate the summation in (|7]) into two parts and change the index of variable as 

LN-l LiV-1 

^2 ^ ^ {LN-k) Cov{wo,L,Wk,L)+ ^ k Cov{wo,L,WLN-k,L)\ (8) 

^ ^ k=l k=l ' 
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Using Q and the fact that the covariance Cov{wo^l, wu^l) in (|^ is symmetric and periodic with 
the period LN, is rewritten as 



684 

Tn 


+ 


LN-l 
2 - ^ ^ 

— ^ Cov(w;o,L,'«;fc,L). 

k=l 


Consider the following two extreme cases using Q. 
1) For the Nyquist Sampling Rate: We have 


a 


2 

1 


684 


2) For the Continuous OFDM Signal: We have 


LN-l 

2 


Cov{wo,L, Wkp) 
72 


k=l 


LN-l 

2 


AT iini X 

IStI L-s-oo 

k=l 


72 

36 


N TT 
2 


9 


sin(a;) 

NM%) 

117 2799 


sin(x) 

2 

+ 9 


+ 9 


sm X, 


iVsin(f) 


sin(^) 

iv^h^ 

4 

+ 


+ 


sm X 


iVsin(f) 


sin(^) 


dx 


SATS JV3 


SAT 


TT 

z 


(9) 


Unless is too small, it becomes approximately cr^ ~ 2799/577. The detailed derivation of 
the integration in Q is explained in Appendix. 


B. Distribution of 

1) The Nyquist Sampling Rate Case (L = 1): In this case, tCn.i’s are Ltd. and thus is 
asymptotically Gaussian distributed due to CUT. That is, 

where a.s. means that the random variable is asymptotically distributed. 

2) Natural Number Sampling Rate Case (L is a natural number): It is easy to see that in 
this case, is asymptotically Gaussian distributed. We divide LN OFDM signal components 
Snp into L subsets 5o, 5i, ■ ■ ■ , Sl^i in the interleaved pattern as 

Sq = {sn,L \ n = Lp + q, 0<p<A^-l}, 0<g<L-l. 
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Since the eomponents Sn,L in the set Sg are mutually independent, CLT ean be applied to eaeh set. 
Then we have L asymptotieally Gaussian distributions that are eorrelated. The sum of eorrelated 
Gaussian random variables is still Gaussian. Therefore, when L is a natural number, we have 



3) Continuous OFDM Signal Case (L —)■ ooj.‘ Clearly, the distribution of ^oo eonverges to 
some distribution as L inereases and we already eheeked that the distributions of for all 
natural number L are asymptotieally Gaussian distribution. Therefore, the distribution of ^oo is 
also Gaussian sueh that 



4) L is a Real Number: Under some assumptions, the above three eases ean be integrated 
into one general result. The CLT ean be applied to m-dependent random proeess [UHl, where 
m-dependent means that two samples from a random proeess with the interval larger than m 
have no statistieal dependeney. Wn,L is an m-dependent random proeess from the faet that the 
eorrelation of the stationary random proeess Wn,L may rapidly diminish as the interval exeeeds 
Ts/N [|71. Therefore, for any real number L larger than one, S,l can be eonsidered as Gaussian 
distributed. 

C. Distribution o/RCM^I^g 

For the Gaussian distributed S,l with mean 6 and varianee eomplementary CDF (CCDF) 
of is given as 



Thus, CCDF of RCM^ L^ is 

^ idB 


F(RCML|^g > a) = P{^L > 10 “o 



where erf (a) = ^ Jq ^ dt is the error funetion. 
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Gaussian approx. (A^=512) 
Gaussian approx. (A^=1024) 
Gaussian approx. (A^=2048) 
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Simulation (A^=512) 
Simulation (A^=1024) 
Simulation (A^=2048) 
Simulation (A^=4096) _ 



dB 


Fig. 1. Comparison of simulated and analytical CCDFs of RCMoo 


Fig. [^compares the simulated and analytieal CCDFs of RCMoo Though we have obtained 
the distribution of RCM^^j^g for general value of L, due to laek of spaee, we only present the 
comparison when L —)■ oo, which is of importance practically. In the simulation, L is set to 
32, which is enough to represent the continuous OFDM signal, and 16-quadrature amplitude 
modulation is used. It is widely known that the metrics describing the envelope behavior of 
OFDM signals such as CM and PAPR do not depend on the modulation order. Note that analysis 
is based on Gaussian approximation in Section IV-B[ Unfortunately, Weibull distribution is a 
kind of heavy-tailed distribution and thus the sum of Weibull random variables slowly converges 
to Gaussian distribution as N increases but it shows good agreement when N is large as shown 
in Fig. 
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V. Oversampling Rate to Capture the CM of Continuous OFDM Signal 

Radically, the metrics such as PAPR and CM have to be calculated from continuous OFDM 
signals. However, it is impossible to handle with the continuous OFDM signal and instead one 
calculates the metrics from the oversampled discrete OFDM signal with sufficient oversampling 
rate. In this case, the metrics calculated from the oversampled OFDM signal can be viewed as 
an estimator for the metrics of the continuous OFDM signal. For instance, in the case of the CM 
calculation, RCM^ is the estimator for true parameter RCMoo. It is natural that a high sampling 
rate guarantees low estimation error but it also entails high complexity. Thus, finding a sufficient 
sampling rate is of great importance. 

In [IHl, upper bounds on the estimation error of PAPR according to L is derived and it is 
proposed that four times oversampling rate is enough to capture the PAPR of continuous OFDM 
signals. However, in the case of CM, this approach is not useful. Differing from the PAPR, the 
estimator RCM^ rarely shows extremely large estimation error. Therefore, the total inspection 
approach is not suitable in the case of CM. This is due to the fact that the Weibull distribution 
is a kind of heavy-tailed distribution which means that its variance is quite large. For instance, 
when the OFDM signal sequence is an impulse signal with N = 1024, it has RCMi = 1024 
and RCMoo — 759. In this case, the square error |RCMi — RCMooP ~ 7 x 10^ is much larger 
than the mean square error (MSB) E[|RCMi — RCMooP] — 5 x 10“^. Thus, in this section we 
statistically approach this problem by deriving the MSB E[|RCMi — RCMooP] according to L. 


A. Joint PDF of 

To obtain the MSB E[|RCMi 


RCMo 


= E[|a/|l — a/^P], first we find the joint PDB 


of and ^oo which is clearly bivariate Gaussian distribution from our investigation in Section 


IV-B We already checked their mean values E[,^l] = E[,^oo] = 6. Next, the correlation coefficient 
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between and is given as 
Cov(^L,^oo) 


PiL,ioa 


^L^OO 


= lim 

L'^OO 


= lim 

L'^OO 


\~^LN—1 \~^L' N—1 ( \ 

Ln=0 Ln'=0 Cov(Wn,L,Wn',LO 

aLCT^LUm 
E'fo"' Ei=o“' 36(9p? + V, + 


LN L'N 


aLCT^LL'N^ 


36 


LN-l ^_i^+Ar^ 


aidooT^LN'^ 


E 

n=0 


nrr 

L 


sinx 


Afsm(f) 


+ 9 


sinx 


Nsm(§) 


sinx 


NM§) 


dx 


where the equation in the integration is periodic with the period Nn. Thus, we have 


P^LiSoo 


36 




9 


sinl^xj 

Nsm{§) 


+ 9 


sin^xj 

Nsm{§) 


+ 


sini^xj 

Nsm{§) 


dx 


( 10 ) 


aLdoo^N Jq 

_ ^OO 

(^L * 

In terms of estimation theory, = cTooIcl implies that can be considered as an unbiased 

minimum MSB estimator of ^oo- That is, orthogonality principle ]E[(,^i — ,^oo) ■Coo] = 0 and 
unbiased property E,[^l\ = lE[,^oo] are satisfied. 


B. MSB between and 

Using the orthogonality principle, the MSB between and ,^oo is obtained as 


E[\^L-U"] = al 




( 11 ) 
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Analytical {N= 
Analytical (N= 
Analytical (N= 
Analytical (N= 
Simulation (N- 
Simulation {N- 
Simulation {N- 
Simulation (N- 



256) 

512) 

1024) 

2048) 

256) 

512) 

1024) 

2048) 


1.0 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2.0 

Oversampling rate (L) 


Fig. 2. Comparison of simulated and analytical — Coo|^] for various L. 


Fig. [^compares — ^ooP] values obtained by simulation and analysis, whieh shows a 

good agreement. Analytieal results are given by using ( [TT] ) and simulation results are obtained 
by testing randomly generated 10® OFDM signal sequenees. 

C. MSE between RCM^ and RCMqo 

The MSE between RCM^ and RCMoo ean be obtained as 

E[|RCMi - RCMoop] = E[|^ - 

= 12 - 2 ■ E[V^] 

AOO nOQ _ 

12 - 2 / / (ei, 6) di, d^2 

Jo Jo 

where (^i, ^ 2 ) is the joint Gaussian PDF of .^1 and ^2 with the eorrelation eoeffieient 
in ([Tg). 

Fig. eompares E[|RCM 2 , — RCMooP] values obtained by simulation and analysis using 
Gaussian approximation. From a praetieal viewpoint, MSE is usually normalized as E[|RCMl — 
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RCMooP]/IE[|RCMoop] = E[|RCM 2 , — RCMooP]/6. Thus, one can conclude that 1.7 times 
oversampling gives the normalized MSB smaller than about 10“^ for practical value of > 256, 
which is remarkable because the FFT size is typically 1.7 times larger than the nominal bandwidth 
for the LTE cellular communication systems [[191 . 


Analytical (N- 
Analytical (N- 
Analytical (N- 
Analytical (N- 
Simulation (N- 
Simulation (N-- 
Simulation (N-- 
Simulation (N-- 



256) 

512) 

1024) 

:2048) 

256) 

512) 

1024) 

2048) 


1.0 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2.0 

Oversampling rate (L) 


Fig. 3. Comparison of simulated and analytical E[|RCMl — RCMcx>P] for various L using Gaussian approximation. 


VI. Conclusion 

In this paper, the properties of CM are investigated. First, asymptotic distribution of the CM is 
obtained, which shows a good agreement to simulation results when the number of subcarriers is 
large. Second, the oversampling rate good enough to capture the CM of the continuous OFDM 
signal is investigated. We confirmed that 1.7 times oversampling rate is good enough for that 
purpose from a viewpoint of MSE, which is of great importance because the FFT size is typically 
1.7 times larger than the nominal bandwidth for the FTE cellular communication systems. 
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Derivation of 


Nn . / \ 


0 ViVsin(f). 


Appendix 
’’dx When r =2, 4, and 6 


1\ TT 

It is known that there are several methods to derive Jq^ (■ 


sm(x) 


"dx when r = 2. By 


/O VArsin(f). 

expanding one of them, both oases of r = 4 and 6 are derived analogously as follows. Without 
loss of generality, by using N = 2M and ohange of variables, we have more simplified expression 
as 


Ntt 

2 


sini^xj 


dx = 


2^M 


r—1 


sin(2Ma;) 

sin(x) 


dx 


where 


sin(2Ma;) 


smuc 


= 2'’(cos(a;) + cos(3x) + • • • + cos((2M — l)a;))^. 


( 12 ) 


Since the integral is over [0, tt], only part of coefficients contribute to the integration. For example. 


when r = 2, all the cross terms from the polynomial expansion in (12) becomes zero after 
integration, i.e., cos(x) cos(3x)dx = 0. After some manipulations, we have 


/ sin(2Ma;) 

\ sin(a;) 

/ sin(2Ma;) 

\ sin(a;) 

sin(2Ma;) ^ ® 
sinfa;) 


= 2^ — + Z2(x 


= 2^( — (m + 8M3) +z4x) 


= 2®( — ( M + + 44M^ ) + zeix) 


where Zr(x) denotes all terms which become zero after integration over [0,7r]. Finally, we have 

2 


Ntt 

2 


Ntt 

2 


N TT 
2 


smuc 


Nsm(§) 


smuc 


Nsm(§) 


smuc 


7Vsin(f) 


dx=- 

2 


, 1 
dx = 7r\ - + 


6N^ 
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